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The  flow  of  fluid  around  obstacles  in  two  dimensions  is 


described  by  the  compressible  Navier  Stokes  equations 


where 


!H  +  II  +  H  =  o 
at  3x  ay 


U  =r  p 


E  =  rpu 


pu  -  o, 


puv  -  c 


■  pue  -  uo  -  vc  -  q 
xx  xy  x 


vo  -  u c  -  q  - 
yy  '’xy 


oxx  =  -  p  -  (2/3)  m$.u  +  2mux 

c  -  m ( u  +v  ) 
vxy  '  y  x 

o  =  -  p  -  (2/3)  m$.u  +  2mu 

y  y  x 


p  »  PRT  m  =  m(T ) 


q  =  kT 
^x  x 


q  =  kT 

y  y  2  2 

e  =  cvT  +  - — ~~ 


where  the  four  variables  p,  u,  v  e  represent  the  physical 
quantities  of  density,  x-  and  y-  components  of  velocity,  and 


internal  energy. 
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This  nonlinear  system  of  mixed  parabolic-hyperbolic  type  in 
two  space  dimensions  and  time  with  four  independent  variables 
must  be  solved  in  an  exterior  region  in  R  .  The  geometry  will 
depend  on  the  particular  physical  situation  that  one  is 
attempting  to  model. 

The  situation  we  shall  be  interested  in  occurs  in  modelling 
flight  conditions,  in  which  the  conditions  at  infinity  are  pre¬ 
scribed  with  a  large  u-velocity  and  v-velocity  zero  and  pre- 
scibed  pB  and  e^.  Fluid  flows  around  an  obstacle  in  x-y  space 

Usually,  this  equation  is  solved  numerically  using  a  finite 
difference  scheme  of  the  Lax-Weudroff  type,  such  as  the 
MacCormack  ADE  mettrod.  Since  these  calculations  can  only  be 
made  on  a  finite  grid  of  points  in  xy  space,  an  artificial  far- 
field  boundary  is  created.  This  boundary  ought  be  sufficiently 
far  away  from. the  object  around  which  the  fluid  is  flowing  so 
.that  local  phenomena  are  not  omitted  by  omitting  part  of  the 
region  of  fluid  flow.  On  the  other  hand,  the  farther  away  the 
region  is,  the  more  grid-points  need  to  be  included  and  thus  the 
more  expensive  time-consuming  the  computations  become. 

One  then  has  the  problem  of  deciding  what  effect  this  new 
boundary  has  on  the  solution  of  the  problem.  Because  of  the 
viscosity  terms  in  (1)  and  the  additional  artificial  viscosity 
introduced  by  the  finite  difference  schemes,  some  boundary  condi¬ 
tions  must  be  imposed. 

As  we  shall  show  in  this  report  on  numerical  experiments, 
cosiderable  care  must  be  exercised  in  the  choice  of  the  boundary 


cond i t ions  . 


If  one  is  interested  in  steady  state  flow,  then  one 
starts  off  with  a  bad  approximation,  and  hopes  that  the  errors  in 
the  numerical  solution  propagate  out  of  the  region  as  transitory 
disturbances  in  the  physical  variables.  One  then  expects  to 
converge  to  the  steady  state  flow. 

We  shall  show  in  this  paper  that  the  incorrect  choice  of 
boundary  conditions  can  give  rise  to  some  of  the  following 
phenomena:  i)  reflecting  boundary  conditions,  in  which  the 

disturbances  in  the  physical  variables  represented  by  the  dif¬ 
ference  between  the  steady  state  flow  and  the  initial  conditions 
are  not  allowed  to  exit  through  the  far-field  boundary  but 
instead  continue  echoing  up  and  down  the  grid,  and  giving  rifee  to 
spurious  oscillatory  solutions;  ii)  under-specified  boundary 
conditions,  in  which  the  converged  steady  state  solution  may 
depend  on  the  initial  conditions;  i i i )  overspecified  boundary 
.conditions  in  which  large  errors  are  introduced  before  conver¬ 
gence  takes  place. 

In  this  paper,  we  first  discuss  the  theory  for  simple  linear 
hyperbolic  systems  in  one  space  dimension.  We  then  analyse 
several  computational  experiments  in  the  light  of  this  theory  in 
one  space  dimension.  This  show  up  the  phenomena  which  we  wish  to 
discuss  in  their  simplest  setting.  Finally  we  discuss  the  more 
complicated  setting  of  two  space  dimensions  and  point  out  how  the 
phenomena  of  one  space  dimension  can  be  recognized  in  this  more 
complicated  setting. 


A  REVIEW  OF  THE  LINEAR  CASE 


The  method  used  in  the  calculations  which  are  the  subject 
this  paper  is  the  MacCormack  alternating  direction  explicit 
scheme.  (6)  [15].  This  is  a  multistep  efficient  scheme  which 

reduces  in  the  linear  case  to  the  Lax-Wendroff  scheme.  For  a 
diagonal  N  x  N  matrix  A,  this  scheme  approximates  the  equation 
Ut  +  AUX  =  0  by. 


(2) 


n+1  n  .At.  n  n  . 

u.  =  u.  -  A—  (u.  .-u.  .) 
3  3  Ax'  3+1  3-1 


of 


display 

As  usual, 
represents  the 


.2, At,  .  n  _n,n  , 

♦  A  <«>  (vr2Wi’ 


{ j 1 ^  *  < j 1  represents  the  space  step  and  n 

time  steps.  If  we  are  considering  the  equation 


0t  +  AUx  =  0 

on  the  region  {(x,t),  0<x<l,  t>0}  then  the  analytic  solution  is 
determined  by  the  initial  conditions  and  boundary  conditions  at 
* x  =  0  and  x  =  1.  If  the  first  k  eigenvalues  are  positive  and  the 

remaining  N  -  k  are  negative  than  the  quantities  u^ - u ^  must  be 

prescribed  at  x  ■  0  and  uk+x - uN  must  be  prescribed  at  x  =  1. 

Thus  if  Wj  =  (ux,  u2» - uk,  0  0 - 0)  and  WXI  = 

(0,  0, - 0,  uk  +  1 - uN)  then  for  we  1 1-posedness  ,  the  boundary 

conditions  must  be 

WT  =  f(t)  +  B0WrI  at  x  =  0 

W  =  g ( t)  +  BjW  at  x  =  1 . 

This  gives  a  total  of  N  boundary  conditions.  If  either 
kx(n-k)  matrix  B^  and  the  (n-k)xk  matrix  BQ  are  non-zero,  then 
the  boundaries  are  reflecting ,  that  is  a  wave  in  Wj  travelling 
left  to  right  will  be  reflected  as  a  wave  in  WIX  running  right  to 


4 
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left . 

If  a  boundary  is  supposed  to  be  non-physical,  then  it  should 
not  be  reflecting,  since  the  reflections  would  depend  on  the 
location  of  the  artificial  or  numerical  boundary. 

Now  let  us  consider  the  difference  scheme  of  (2).  If  the 
grid  points  are  given  by  {xjJq  with  x  =  0  and  Xj  =  1,  then  it  is 
clear  from  (2)  that  2N  boundary  conditions  are  required.  Thus  we 
must  prescribe  boundary  conditions  in  such  a  way  as  to  least 
affect  the  closeness  of  the  numerical  to  the  analytic  solution. 

This  situation  has  been  the  subject  of  several  papers.  In 
(3],  Gustaffson  and  Kreiss  point  out  the  danger  of  over¬ 
specification.  By  this  is  meant  that  all  u^  are  specified  at 
both  endpoints.  This  might  be  tempting  because  one  might  argue 

that  if  the  boundary  conditions  - -ik  are  given  constants  at  x 

=  0,  then  eventually  these  same  values  will  be  assumed  by  these 
•variables  at  x  =  1.  However,  in  [4]  it  is  pointed  out  that 
convergence  may  or  may  not  occur,  depending  on  whether  the  number 
of  grid  points  is  odd  or  even. 

A  method  which  works  well  as  pointed  out  in  [3]  is  to  impose 
3u . 

13T  (1'fc)  "  0  1  <  i  <  k. 


3u  . 
l 

3x 


(o,t)  =  0 


or  in  the  numerical  scheme 


u .  ,  =  u .  ,  , 
i ,  J  i  ,  »J—  1 


k+1  <  i  <  n. 


1  <  i  <  k . 


u . ,o  =  u.  .  k  +  1  <  i  <  n . 

1  1  f  1 


This  introduces  small  errors  at  the  outflow  but  these  errors 


do  not  propogate  upstream.  This  is  proved  analytically  by  Parter 
HO]  . 

In  (12] ,  many  different  numerical  boundary  conditions  are 
given.  The  conclusion  is  that  upwind  differencing  at  the  point 


of  outflow 


n+1  n  , . . A t . .  . 

uj  ■  uj  +  il(M),uj-ruj> 


is  most  accurate,  although  it  converges  with  the  same  speed  as 


the  previously  discussed  ux  =  0 . 

The  most  serious  error  which  one  could  make  would  be  to 


prescribe  conditions  at  the  wrong  end.  In  other  words,  since  Uj 
is  right  running,  this  would  involve  prescribing  u^  at  x  =  1  and 
imposing  ulx  =  1  at  x  =  0.  This  would  result  in  convergence  to  a 
steady  state  which  depends  on  the  initial  conditions. 


*  THE  NAVIER-STOKES  EQUATIONS  AND  CHARACTERISTIC  VARIABLES: 


We  now  begin  our  discussion  of  the  equations  of  gas 
dynamics.  We  will  neglect  viscosity  for  the  puposes  of  this 
analysis.  We  will  assume  that  the  flow  is  one-dimensional  and 
subsonic  and  that  the  deviations  from  free-stream  solutions  are 
small.  This  will  allow  us  to  neglect  second  order  terms. 


There  are  many  forms  of  this  equation,  but  the  one  most 
suitable  for  the  present  discussion  is 

2H  +  aJLH  -  o 
at  Aax  u 


where 


and 

M~ 1  =  /  1  0  0 

I  -u/p  1/p  U 

V1— )u2  Jl-Y)u  (y-1) 

Here  we  make  the  key  assumption  that  deviations  from  the 
free  stream  are  going  to  be  sufficiently  small  that  we  can  treat 
the  entries  in  the  matrix  A  as  being  approximately  constant  (at 
least  locally).  Denote  these  frozen  variables  by  O-subscript. 


We  then  make  the  substitution 


-1/c' 


1/p  c 
o  o 

1/p  c 
'  o  o 


and  when  this  is  substituted  into  (6)  we  obtain 
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(5) 


aw  aw 

1  +  u  =  0 


at 


o  ax 


aw  aw 

“at  +  (u+co}  ~aw  “  0 


aw  aw 

+  (u  -cl  -~  =  0 


at 


Jo  o'  ax 


W2  = 


u  + 


p  c 
o  o 


W,  =  -u  +  - T- 

3  p  c 


o  o 


Notice  now  how  this  breaks  down  into  two  separate  cases.  On  the 
one  hand,  if  flow  is  supersonic  then  all  wave  motion  is  in  the 
left  to  right  direction.  In  this  case  all  analytic  bound  / 
condition  sought  be  prescribed  at  the  left  hand  side  and  t  ‘y 
numerical  boundary  conditions  prescribed  at  the  right  han  e. 

Since  the  substitution  (4)  is  equivalent  to 

o  -  Ki  +  |*A)IVV 

u  =  1/2(K2-K3) 
p  =  p0C0/2(K2+K3) 

« it  follows  that  prescribing  all  physical  variables  at  the  inflow 
and  prescribing  3p/3x  =  3u/3x  =  3p/8x  =  0  at  the  outflow  is 
legitimate  in  terms  of  analytical  and  numerical  requirements  in 
the  supersonic  case. 

However,  we  must  now  consider  the  case  of  subsonic  flow.  in 
this  case  the  situation  is  completely  different.  Here,  two  of 
the  variables  W^  and  W2  go  left  to  right  with  velocities  u  and 
u+c  respectively,  whereas  one  of  the  variables  runs  right  to  left 
with  velocity  c0-uQ.  While  the  variables  w2  and  w3  have  no  clear 
physical  significance,  yet  it  is  only  be  considering  these  vari¬ 
ables  that  the  full  wave  structure  of  the  equations  (5)  or  (6) 
can  be  understood.  Thus,  one  would  be  led  to  predict,  for  small 
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deviations  from  free  stream  conditions,  that  the  best  boundary 

conditions  would  be  for  an  interval  (o,L) 

dW 

Wj  ( o ,  t )  =  Kj  (  L,  t )  =  0 

dw 

W2  ( o ,  t )  =  K2  (  L,  t )  =  0  (6) 

dW 

( O ,  t )  =0  W3  (  L,  t )  =  K3 

Note  the  curious  aspect  of  these  boundary  conditions.  In 
order  to  prescribe  the  numerical  values  K ^  and  K2,  we  need  to 
know  accurately  all  three  physical  variables  at  some  distance  to 
the  left.  However,  only  the  two  combinations  Kj  and  K2  are  pre¬ 
scribed.  This  can  be  summarized  by  saying  that  while  we  have 
used  all  three  pieces  of  information  upstream,  we  have  done  so  in 
such  a  way  that  one  degree  of  freedom  remains,  thus  allowing  the 
'waves  in  W,  to  exit  without  reflections. 

On  the  basis  of  the  linearized  model,  various  other  combi¬ 
nations  would  be  well-posed.  For  example,  it  is  possible  to 
prescribe  K3  in  terms  of  either  K-^  or  K2  at  the  outflow  x  =  L. 
Thus  at  the  outflow  one  may  prescribe 

W2  ( L, t )  =  F3  (t)  +  CjW  ( L, t )  +  c2W2  ( L , t ) 

For  example  if  c3  =  0,  c2  =  1,  then  this  amounts  to  putting 
u  ( L, t )  =  1/2  F3  (7) 

i.e.,  we  prescribe  velocity  at  the  outflow. 

Alternatively,  we  might  take  c^  =  0,  c2  =  -1  and  we  would  get 
p  (L,t)  =  ((p0c0)/2)F3  (8) 

i.e.  we  prescribe  pressure  at  the  outflow.  Many  other  combi 


t 
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nations  are  possible,  but  as  remarked  in  section  II,  all  these 
will  cause  errors  in  the  initial  data  to  be  reflected  back  into 
the  medium  as  waves  running  from  right  to  left.  For  example,  we 
would  predict  that  an  error  in  would  be  reflected  back  as  an 
error  in  Wj  if  we  use  boundary  condition  (8).  As  we  shall  see, 
this  is  exactly  what  happens. 

At  the  inflow  end,  we  may  prescribe  and  W2  in  terms  of 

.  Thus  the  following  boundary  conditions  are  well  posed; 

Wj  ( o ,  t )  =  Fj  +  c :  W-j  {  o  ,  t )  (9) 

W2  (o,t)  =  F2  +  c2W3  ( o , t ) 

For  example,  choosing  c2  =  +1  in  (lib)  corresponds  to 

u(o,t)  =  <1/2)F2 

(i.e.  prescribing  u  at  the  inflow)  and  c2  =  -1  corresponds  to 

P(o,t)  =  <P0co/2)F2 

(i.e.  prescribing  p) .  One  can  prescribe  the  combination  (u,p)  by 

first  choosing  c2  =  1  (thereby  prescribing  u)  and  then  choosing 

Cj  =  PQ/cof  thereby  prescribing  p  in  terms  of  a  given  F,  and  a 

prescribed  u(o,t).  Since  (11a)  and  (lib)  reduce  to 

u(l+c„)  +  — j—  (1— c_)  =  F_ 

2  ‘>oco  2  2 

p  -  <Po/Co  -  Cj)  .  CjU  =  Fj 

o  o 

about  the  only  condition  we  cannot  prescribe  is  u(o,t),  p(o,t), 
since  there  is  no  choice  of  c^,  c2  to  eliminate  p  from  these 
egua  t ions . 

Again,  we  emphasize  that  each  of  these  boundary  conditions 
is  reflecting ,  i.e.  deviations  from  the  free  stream  in  the 
initial  data  get  reflected  back  as  waves  in  W-^  and  W2  and  then 


.1 


travel  back  downstream.  About  the  worst  thing  that  can  be  done 
is  to  prescribe  reflecting  boundary  conditions  at  the  inflow 
x  =  0  and  the  outflow  x  =  L.  In  this  case  errors  can  keep  being 
reflected  up  and  down  the  region,  never  being  allowed  to  exit. 
This  prevents  convergence  to  a  steady  state  and  may  even  give 
rise  to  fictitious  periodic  oscillations. 

We  conclude  this  section  with  a  review  of  the  conclusions  on 
boundary  conditions.  Based  on  the  linear  model,  boundary 
conditions  (8)  seem  optimal.  Any  other  prescription  of  the 
physical  variables,  although  well  posed,  casues  reflections  of 
the  deviation  from  the  true  solution.  If  for  example,  only  the 
physical  variable  p^  is  known  at  the  outflow  then  it  is  possible 
to  prescribe  p  at  the  outflow,  in  such  a  manner  that  the  problem 
remains  well  posed.  We  emphasize  that  this  will  cause  errors  to 
propagate  upstream,  thereby  slowing  the  process  of  convergence  to 
.  steady  state.  This  may  not  be  too  bad,  so  long  as  the  upstream 
boundary  conditions  are  not  also  reflecting.  On  the  other  hand, 
if  they  are,  then  convergence  to  free  stream  may  never  occur. 

VI  .  DISCUSSION  OF  NUMERICAL  RESULTS  -  NONLINEAR  COUPLING: 

The  one  dimensional  Na v ier-Stokes  equations  were  solved  with 
an  alternating  direction  explicit  MacCormack  scheme,  on  a  one¬ 
dimensional  net  with  forty  grid  points.  The  code  was  an  exact 
one-dimensional  version  of  a  three-d imens ionsa 1  code  which  had 
proved  successful  in  many  supersonic  studies  [6],  [11].  There 
were  two  questions  to  answer.  The  first  was,  given  that  equation 
(7)  is  correct  for  infinitesimally  small  deviations  from  a 


constant  free  steam,  how  correct  is  it  when  deviations  of  an 
intermediate  (of  the  order  of  10%)  magnitude  are  present 
instead?  It  would  be  too  much  to  hope  that  the  variables  ,  W2, 
and  ,  remain  uncoupled,  but  we  should  be  able  to  get  an  idea  of 
the  order  of  magnitudes  of  the  coupling  involved.  The  second 
problem  of  course,  is  to  assess  the  influence  of  the  various 
types  of  boundary  conditions  commonly  employed.  We  will  deal 
with  the  latter  problem  in  section  VII. 

To  do  this,  we  considered  a  uniform  free  stream  situation, 
with  pressure  equal  to  2000  lbs/ft^,  velocity  equal  to  548  ft/sec 
and  density  equal  to  0.0023  slugs/ft^.  We  created  a  deviation 
from  the  steady  state  condition  in  a  variety  of  ways  as  in  (9), 
flO],  and  then  watched  the  progress  (or  lack  of  it)  to  a  steady 
state.  A  variety  different  wrong  initial  conditions  were  used. 
One  type  was  to  impose  a  10%  deviation  in  one  of  the  variables 
Kj,  K2,  K-j  at  the  points  {x^,  j  =  18 , 19 , 20 , 21 , 22  }  .  We  should  then 
watch  the  disturbance,  graphically  as  it  propagated  up  or  down 
stream.  Another  possiblity  was  to  put  in  uniformly  wrong  initial 
conditions  where  some  of  all  of  the  characteristic  variables  w^ , 
W2 *  w3  are  perturbed  throughout  by  a  percentage  error  of  10%. 

Each  plot  then  showed  the  percentage  error,  with  the  different 
curves  representing  the  progress  of  time  as  one  ascends  the  plot. 
The  curves  are  plotted  every  twenty  five  time  steps 
when  At  =  ( .9) ( Ax)/(u+c) .  We  also  point  out  the  percentage 
1#  W2,  Wj  at  the  end  of  the  run,  so  as  to  obtain 


errors  in  W 


information  on  relative  accuracy  and  speed  of  convergence  of  the 
various  methods. 

Figure  1  gives  the  results  of  an  experiment,  which  is  ideal 
in  terms  of  the  linear  theory.  An  initial  disturbance  in  the 
left-running  characteristic  variable  is  given  and  we  observe 
deflections  in  the  variables  W^,  W2 ,  .  As  one  can  readily  see 

from  the  pictures,  the  disturbance  propagates  upstream  rapidly 
until  convergence  is  reached  (.1%  agreement  with  physical 
variables),  which  takes  place  within  130  time  steps.  This  gives 
us  six  curves.  Note  that  although  a  good  deal  of  undershoot  and 
overshoot  in  W3  becomes  apparent,  there  is  no  significant 
interaction  with  or  W2 .  The  same  situation  appears  with 
initial  disturbances  in  ,  the  slow-moving  right  running  wave. 
From  this  experiment  it  appears  that  deviations  in  either  Wj  or 
will  not  effect  either  of  the  other  two  variables.  However, 
as  shown  in  Figures  2  and  3,  when  initial  disturbances  are  in  W2 , 
a  different  situation  exists.  Figure  2  shows  a  uniform  initial 
disturbance  in  W2  of  minus  ten  percent,  while  and  w3  are  left 
undisturbed.  Initially  the  wave  in  W2,  propagates  rapidly  out  of 
the  medium.  Indeed,  after  fifty  time  steps  it  is  essentially 
gone  from  the  picture.  However,  this  does  no  happen  without 
affecting  the  of  the  two  variables.  Notice  how,  in  the  top 
graph,  a  large  disturbance  is  left  in  after  fifty  time  steps 
and  in  W3,  we  have  that  W3  values  one  almost  constant  at  minus 
eight  percent.  However,  once  the  W2  wave  has  made  its  exit,  the 
other  two  variables  uncouple  and  resume  their  normal  wave  motion, 
and  the  error  can  be  seen  propagating  out  of  the  solution  in  the 


usual  way  convergence  is  attained  within  three  hundred 
interations.  These  particular  results  illustrate  what  became 
increasingly  clear  throughout  the  series;  perturbations  in  W2  had 
a  large  effect  on  and  whereas  if  W2  was  not  perturbed, 
and  W-j  behaved  as  if  uncoupled.  Perturbations  in  and  had 
little  effect  on  W2 .  No  gualitative  explanation  of  the 
phenomenon  is  know  at  this  time.  Figure  3  shows  the  same 
phenomenon.  Here  an  error  of  -10%  is  made  in  W2  whereas  errors 
of  +10%  are  made  in  Wj  and  W2 .  Again,  we  see  that  until  the  W2 
wave  exits,  there  are  massive  disturbances  in  the  wave  structive 
of  Wj  and  W2 .  As  soon  as  W2  exits,  (after  75  interations)  the 
regular  wave  structive  reasserts  itself  and  errors  propagate  out 
in  predictable  wave¬ 
like  manner.  Again,  convergence  takes  approximately  three 
hundred  and  -twenty  five  iterations.  It  seems  clear  that  this  is 
optimal  given  the  limited  wave  velocity,  so  we  can  deduce  that 
these  affects  are  due  to  the  nonlinear  coupling.  Thus,  even 
after  the  W2  wave  exists  it  will  take  at  least  a  maximum  time 
of  {L/(u-c),  L/u }  seconds  for  the  resulting  errors  to  propagate 
out  of  the  system. 
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VII.  DISCUSSION  OF  NUMERICAL  RESULTS  -  BOUNDARY  CONDITIONS: 

The  linear  "small  deflection"  theory  predicts  that  the  best 
boundary  conditions  would  be  the  prescription  of  the  character¬ 
istic  variables  at  their  point  of  entry  with  some  form  of 
(stable)  numerical  boundary  condition  for  the  point  of  exit. 

Here  are  two  such  schemes 

INFLOW  OUTFLOW 

dW. 

Wj  ( O  r  t )  =  K1  (L,t)  =  0 

dw 

W2  (o,t)  =  K2  (  L,  t )  =  0  (12) 


dw 

(0,t)  =0  W3  (  L,  t )  =  K3 

Here  and  K2  are  numbers  calculated  from  the  known  values 
of  u,  p,  p  at  the  inflow  and  K-j  is  calculated  from  the  known 
values  of  u  and  p  at  the  outflow.  Notice,  however,  the  one 
"degree  of  freedom"  is  left  at  the  inflow  point.  This  allows  the 
variables  to  adjust  but  in  compensating  ways.  The  boundary 
conditions  in  the  code  are  usually  in  terms  of  the  physical 
variables  so  we  translate  (12)  to  physical  variables. 

INFLOW 

iK2  U2  -  (1/‘,oCo>P21 

Uj  =  1/2  [Kj  .  u2  -  (l/p0c0)p2J 

P1  =  Ki  +  (po/2co)[K2  "  U2  +  {1/poco)p21 

OUTFLOW 

UN  *  1/2|un-l  +  (1/»oco)  PN-1  +  K3 

PN  ’  <VV2)|K3  +  UN-1  +  (V»oeo'  PN-1] 

CN  *  ,l>c/2co,^K3  *  uN-l'  ■  [1/(2cq)  ]  PN_J  +  PN_J 
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This  set  of  boundary  conditions  is  predicted  to  work  well  in  the 
linear  studies  of  one  equation,  occurring  in  (21  and  (31 .  We 
shall  call  these  boundary  conditions  the  "no-change  character¬ 
istic  boundary  conditions".  Another  possiblity,  suggested  by 
one-D  analogues  in  (21,  is  the  following: 


INFLOW 


W2  ( o ,  t )  =  Kj 


w2  ( O , t )  =  k2 


OUTFLOW 

3W, 


IT  +  UN 


3W1 
3  x 


3W2 
~ jft 


=  0 


3  W, 


+  (u-c) 


N  3  x 


=  0  (13) 


3W3 

3t 


+  (u-c) 


3W3 
1  3x 


=  0 


W3  ( L,t)  =  K3 


where  derivaties  in  the  x-variable  are  downwind  at  the  inflow  and 
upwind  at  the  outflow  and  forward  in  time.  The  numbers  Kj ,  k2, 

K3  are  prescribed  as  before.  In  terms  of  the  physical  variables 
these  translate  into 


INFLOW 


u 


n  +  1 
1 

n  +  1 


=  1/2 [K,  + 


U1  •  11/poco)pl  * 


/  _  n  n 

(uo'co)(i7>  luru2  + 


2  ”1  "o''o'*'l  ’  '-o  ''o,'Ax,,v'l  “2  ^1 

( poco/2) fK2+( 1/poco)pl-ul{uo"co) ( At/Ax)  [u2ul+  ( 1/poco) (Pi+P2J 1 1 


p?)  n 


pl  +  1=  Kl+(po/2co} rK2+(1/poCo)pl"Ul+(uo"co) (At/Ax) [u2-ul  +  (1/poco) (pi"p2)1  ’ 


OUTFLOW 


=  (1/2)[UHVT  PN"K3  +(At/Ax)(uo+Co)  [UN-rUN+(1/poCo)(pN-l  +  PN)n 
o  o 


n 

V 

n 


n  1 
)  J  1 

n  n, 

n, 


PN  =  '»oVJ)  |K3+UN+<1/‘,oC0)PN+(4t/4x)(Uo+Co)  IUN-rUN*<1/t'oCol(PN-rPNl 
PN+1=  (Po/2CO,|K3+UN+<1/(,oCO)PN+(4t/4X),UO+Co)UN-rUN+1/0OCOPN-rPNn 
+  CK  '  U/co,pN  +  (4t/4x'uolpS-l-p13+U/c2)(P^  -  Pk_j >  1 
We  shall  call  these  the  "windward  difference  characteristic 

boundary  conditions".  Note:  u0,  cQ  can  be  different  values  at 
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inflow  and  outflow.  The  performance  of  the  code  with  either  of 
these  was  analyzed  by  posing  initial  conditions  in  which  there 
was  a  disturbance  in  one  or  more  of  the  characteristic  variables 
either  locally  at  the  center  of  the  grid  or  uniformly  throughout 
the  grid,  of  the  order  of  10%. 

Thus  figure  1  shows  what  happens  if  the  disturbance  in  only 
in  the  third  characteristics  variable  locally  using  the  windward 
characteristic  variables. 

Figure  3  shows  the  effect  of  a  plus  +10%  error  in  the 
initial  conditions  Wj  and  W3  and  a  -10%  error  in  W2 .  (The  first 
curve  from  the  bottom  is  the  initial  state  of  th  variable,  and 
the  others  are  the  states  at  intervals  of  25  intera t ions )  .  In 
figure  2,  we  have  an  initial  disturbance  in  W2  of  -10%  with  no 
initial  disturbance  in  or  W3 .  The  pictures  look  essentiallly 
the  same  as  figure  1.  There  is  considerable  nonlinear 
interaction  until  the  W2  wave  exists,  and  then  uncoupled  wave 
motion  to  the  right  in  the  first  variable  (W3)  and  to  the  left  in 
the  third  variable  (W3).  There  are  no  reflections  when  the 
and  W3  waves  exit  and  convergence  is  reached  in  300  iterations. 

These  computations  were  made  using  the  windward  differencing 
characteristic  boundary  conditions,  although  the  same  results 
were  obtained  with  the  no  change  characteristic  boundary 
conditions . 

In  Figure  4,  we  show  the  effect  of  a  local  disturbance  at 
the  center  of  the  grid  in  the  W2  variable  with  the  second  set  of 
B.C.'s  and  in  figure  5  we  show  a  speeded  up  version  (every  50 
iterations)  of  the  same  disturbance  with  the  first  set  of  B.C.'s. 


The  third  and  fifth  graphs  on  figure  4  are  almost  exactly  the 
same  as  the  second  and  fourth  on  figure  5.  Figure  5  shows  con¬ 
vergence  being  reached  in  300  iterations. 

We  conclude  that  either  of  the  first  two  sets  of  boundary 
conditions  give  optimal  convergence  since  convergence  cannot  take 
place  until  the  wave  in  W2  exists  (very  quickly)  and  the  residual 
(nonlinear)  effects  of  Wj  on  W-j  can  exit  upstream.  If  they  can 
do  this  without  any  reflections,  then  the  convergence  is  essen¬ 
tially  optimal . 

Sometimes,  it  is  objected  that  in  a  wind  tunnel  experiment, 
the  only  variable  known  downstream  is  pressure  and  that  we  are 
requiring  too  much  information  in  prescribing  ,  which  demands  a 
knowledge  of  p  and  u  at  the  outflow.  Suppose,  then,  we  just  pre¬ 
scribe  pa  at  the  outflow  using  the  otherwise  successful 
3W2  3W2 

conditions  of  — —  =  =  0  as  complementary  numerical  boundary 

*  conditions.  Then  we  could  have,  for  example 


INFLOW 


Wj  ( o , t )  =  Kj 


W2(o,t)  =  K2 


OUTFLOW 

P  =  P„ 


dw. 


6Vll 

dx 


dw. 


=  ( o ,  t )  =  0 


(  L, t )  =  0 


=  (  L, t )  =  0 


(  14) 


dx  -  w  ~dx 

The  inflow  boundary  conditions  are  precisely  those  of  (12). 

outflow  boundary  conditions  (used  by  Steger  [12])  are 


The 
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PN  =  P. 

PN  -  ll/cj]|p.  -  PN_! )  +  PH_! 

UN  '  (1/coco)(Pn-1  ‘  p->  +  UN— 1 

The  predictions  of  section  VI  are  clear.  The  fact  that  p  is 
prescribed  means  that  when  a  wave  in  W2  comes  downstream,  it 
exits  by  adjusting  the  u  values  a  c  x  =  xN.  This  in  turn  causes 
disturbances  in  W_=  -u  +  (1/p  c  )p  which  cause  a  reflected  wave 
upstream.  This  wave  can  be  seen  by  comparing  figure  (6)  with 
figure  (4).  Since  the  upstream  B.C.  is  non-reflecting,  this 
means  that  the  left  running  wave  will  exit  without  incidence. 

When  composed  with  boundary  conditions  (12)  or  (13)  it  is 
obviously  less  desirable  because  of  the  magnitude  of  the 
reflection  in  W3 .  However  it  does  converge  in  approximately  300 
iterations,  which  is  again  almost  optimal.  The  effect  of  these 
large  oscillations  in  more  complicated  geometries  may  prove 
undesirable,  however. 

We  briefly  review  our  progress  so  far.  Two  sets  of  non¬ 
reflecting  boundary  conditions  have  been  produced.  Doth  of  which 
give  optimal  convergence  but  which  rely  on  a  great  deal  of  infor¬ 
mation  at  both  ends.  The  information  however  is  used  in  such  a 
way  as  to  allow  additional  degrees  of  freedom  for  the  waves  to 
exit  without  repeated  reflections.  One  reflecting  and  non¬ 
reflecting  boundary  condition  can  be  combined  to  obtain  almost 
optimal  convergence  at  the  cost  of  some  large  left  running 
reflections,  whose  effect  in  more  complicated  geometries  remains 
uncertain.  While  B.C.  (12)  was  expected  to  be  less  accurate  than 
B.C.  (13),  little  evidence  for  this  has  been  uncovered,  except  at 
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the  boundaries.  Time  dependent  periodic  flows  (e.g.  self  excited 
oscillations)  may  prefer  B.C.  (13)  however.  We  now  consider  some 
of  the  other  boundary  conditions  which  have  been  tried  previously 
in  the  literature. 

First  we  consider  the  case  of  reflecting  boundary 
conditions .  These  occur  in  several  places  in  the  literature,  for 
example  in  [10]  and  [12]  and  have  been  discussed  in  section  VI. 

As  we  have  seen,  these  arise  from  prescibing  combinations  of 
characteristic  variables  such  as  pressure  downstream,  and  other 
combinations  (perhaps  to  density  and  velocity  upstream).  In  [10] 
Rudy  and  Strikwerda  considered  (among  many  others)  the  boundary 
conditions 


INFLOW 
u  =  u 

QO 

T  =  T 

00 


and  in  [12]  Steger  uses 


OUTFLOW 


du 

dx 

dp 

dx 


0 

0 


P  =  Pco 


dw1 

dx 


0 


(15) 


dw 


p  =  »„ 

dw 

-di  =  0 


i 


dx 
P  =  P 


=  0 


(16) 


Figure  7  show  the  effects  of  an  initial  local  disturbance  in  W-^ 
on  both  of  these  sets  of  boundary  conditions,  (13)  on  the  left 
and  (14)  on  the  right.  Note  that  first  there  is  only  a  distur¬ 
bance  in  the  bottom  picture.  By  iteration  75  (fourth  curve  up 
from  the  bottom)  we  can  see  reflections  in  both  W2  and  W^ 
although  the  W2  deiration  is  a  little  harder  to  see.  By 
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iteration  125  it  can  be  seen  that  the  W2  wave  has  travelled  down¬ 
stream  and  is  reflected  back  upstream  in  .  These  reflections 
continued  (somewhere  smeared  out)  for  a  least  2000  iterations. 
Perhaps  most  startling  is  (at  least  at  the  beginning)  how  similar 
they  are.  A  conclusion  may  be  drawn  from  B.C.'s  (14),  (15),  and 
(16).  Prescribing  physical  values  instead  of  characteristic 
values  gives  rise  to  reflections.  If  the  reflections  can  occur 
only  at  one  end,  this  does  not  impede  convergence.  However, 
reflecting  conditions  at  both  ends  can  be  disastrous.  This 
accounts  for  the  "spurious  pressure  waves"  mentioned  by  Moretti 
17]  . 


The  case  of  prescription  at  the  wrong  end  is  now  considered. 
In  this  case,  the  variables  u,  p,  p,  are  prescribed  at  the  inflow 

and  the  conditious 

du  _  d£  _  djp  _  n 
dx  dx  dx 

are  prescribed  at  the  outflow.  In  this  case,  square  wave  distur¬ 
bances  which  effected  only  the  interiors  of  the  domain  exited  as 
in  figure  1,  with  some  minor  oscillations.  However,  when  a  uni¬ 
formly  wrong  initial  condition  was  imposed,  convergence  was  very 
slow  with  large  oscillations,  and  the  solution  converged  to  the 
wrong  values,  with  errors  of  as  much  as  27%.  The  converged  value 
was  a  function  of  the  initial  condition  at  the  outflow  end,  as 
predicted  in  Gustafson  and  Kreiss  (31 .  The  resulting  graphs  are 
given  in  figure  8. 

Related  to  the  above  problem  is  the  method  of  over 
prescription  of  boundaries.  This  method  is  mentioned  in  [9]  as 
giving  good  results  although  the  authors  caution  against  it  on 


the  grounds  of  small  oscillations  being  present.  In  fact  the 
situation  is  much  more  serious.  If  initial  waves  in  the  interior 
of  the  domain  are  used  with  the  initial  conditions  correct  near 
the  outflow,  then  the  solution  converges  rapidly  as  the  travel¬ 
ling  waves  exit  without  reflections  and  with  minor  oscillations. 
However,  if  the  initial  data  is  uniformaly  wrong  with  an  error 
initially  at  the  outflow  point,  the  and  Wj  converge  rapidly 
but  accumulates  huge  errors  of  the  order  of  70%.  Eventually 
when  Wj  and  W2  are  converged,  the  correct  value  is  propagated 
upwind  in  but  taking  large  amounts  of  time  to  converge  because 
of  the  large  errors  near  the  inflow  point.  Indeed,  as  W2  becomes 
more  accurate  downstream  the  inaccuracies  become  much  larger 
(140%)  upstream.  Particularly  in  a  time  dependent  problem  or  in 
a  problem  with  more  complicated  geometries,  this  could  be  truly 
disastrous.  It  points  to  another  fact:  if  a  new  boundary 
‘condition  is  being  tested,  it  is  not  sufficient  to  consider 
initial  value  perturbations  from  free  stream  which  are  non-zero 
in  the  interior  only.  In  this  case,  we  might  have  drawn  totally 
wrong  conclusions  from  the  time  taken  for  convergence. 

In  [10],  a  separate  non-reflecting  boundary  conditrion  is 
proposed.  This  boundary  conditon  alters  the  value  of  K^, 
increasing  it  if  the  computed  value  of  p  is  less  than  pB,  and 
decreasing  it  if  the  computed  value  of  p  is  greater  than  p^ . 

Some  thought  shows  that  this  cannot  be  optimal.  Indeed,  we  can 
choose  a  variation  from  the  initial  conditional  in  which  p  is 
less  than  p^,  but  because  u  is  smaller  than  ub,  the  computed 
is  actually  larger  than  W_  .  In  this  case,  the  proposed  non- 


FIGURE  4 

A  square  wave  in  W2  exits,  leaving  trailing  disturbances 
but  no  reflection.  See  figure  (6)  if  outflow  B.C.  is  reflecting 
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TWO  DIMENSIONAL  RESULTS 

Having  now  understood  the  phenomena  which  can  occur  when 
calculations  are  made  with  one  space  dimension,  we  now  consider 
the  vastly  more  complicated  situation  of  two  space  dimensions. 

In  this  situation,  we  shall  solve  the  Navier  Stokes  equation 
(1),  again  by  the  standard  MacCormack  A.D.E.  scheme  on  a  20  x  20 
grid,  which  is  a  very  simplified  model  of  a  wind  tunnel.  We 
shall  continue  to  consider  flow  close  to  the  free  stream  flow 
previously  studied  in  the  one  dimensional  case.  The  physical 
values  are  given  earlier.  Figure  8  shows  the  geometry  of  the 
situation.  The  fluid  is  flowing  in  at  the  top  right  corner  of 
the  grid  and  flowing  out  at  the  bottom  left. 

We  have  three  distinct  types  of  boundaries  to  consider.  We 
have  the  inflow  and  outflow  boundaries  (as  before)  and  in 
addition,  two  sidewall  boundaries,  where  the  fluid  is  flowing 
parallel  to  the  boundary. 

We  shall  continue  to  impose  one  dimensional  boundary 
conditions  of  the  type  given  in  the  first  section  on  the  inflow 
and  outflow,  along  with  the  additional  condition  v  =  0.  This 
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where 


A  =  r0 


~2"(  Y~1 )  -u(y-1)  0 


|y-^(y-D  o 


Freezing  the  coefficients  of  A  in  (6),  we  substitue 

t  =  p  +  — £— 

1  2 
CO 


T  =  v  +  — p 

2  p  c  * 

Ko  o 


T.  =  -v  +  — - — p 
3  »oco 

T4  *  u 

Equation  (6)  is  then  transformed  to 
Tit  -  0 

<8>  T2t+CoT2y=° 

T3t  -  co  T3y  =  0 
T4 1  =  0 

Thus,  if  we  only  consider  deviations  transverse  to  the  free 
stream  flow,  we  have  on  the  basis  of  the  linearized  model,  four 
non-physical  variable  T^,  T2 ,  T3,  T^,  two  of  which  u  and  entrophy 
Tj  move  with  zero  velocity  in  the  y  direction,  one  of  which,  T2, 
moves  with  speed  cQ  in  the  positive  x  direction,  and  one  of  which 
moves  with  speed  cQ  in  the  negative  x  direction.  Thus  one 
dimensional  theory  predicts  that  at  the  sidewall,  we  should 
impose 


•  A 
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Top  wall 


P  c 
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[v  .  .  +  — ■ * - 

J"1  poco 


l/2k3  +  l/2(v._1  +  Pj_1/P0cQ) 


’j  -  T^T  +  ^  - 

J  o  J  Ko  o 

uj  =  uj-l* 


k3>  +  "j-1  -  ~2  Pj-: 


Bottom  wall 


U1  =  u2 

V1  =  !/2 1  v2  +  k2  “  P2/(pocon 

P  C 

P1  =  '°T  {“V2  +  k2  +  P2/(poCo)J 

P1  =  p2  +  (po/2Co)  [_V2  +  k2  "  P2/poCo] 
where  the  free  stream  values  are  used  to  specify  the 


characteristic  combinations 

Poo 

w2  =  v- +  7T  k3  =  v„,  -  pyp0c0- 

o 

DISCUSSION  OF  TWO-DIMENSIONAL  COMPUTATIONS 


We  now  discuss  a  series  of  computations  which  illustrate  the 
applications  of  the  boundary  conditions  we  have  developed,  as 
well  as  their  short  comings.  To  illustrate  varying  boundary 
conditions,  we  shall  consider  the  experiment  where  initial 
conditions  are  imposed  to  create  a  10%  error  in  T2,  (the  left¬ 
ward  moving  variable),  whereas  Tj_,  T^,  T^  are  at  the  correct  free 
stream  constants.  The  ten  percent  error  is  imposed  for  all  x  at 
five  central  y  grid  values.  Thus,  if  the  linear  model  was 
exactly  accurate,  the  square  wave  would  exit  in  the  positive  y- 


direction  with  velocity  cQ. 
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In  figure  9,  this  is  shown  happening.  Figure  9  has  Neuman 
boundary  conditions  imposed  at  the  inflow  and  outflow,  to 
illustrate  the  case  of  flow  in  an  infinite  region.  Thus,  there 
is  no  contribution  to  the  deviation  from  free-stream  from  the 
inflow  and  outflow.  The  picture  looks  almost  identical  to  the 
one-dimensional  plots  where  the  y-variable  was  ignored.  There  is 
some  overshoot  visible  but  the  wave  motion  is  as  predicted  and 
the  disturbance  in  the  variable  u  +  — —  p  exits  uneventfully, 

DC 

o  o 

leaving  some  trailing  disturbances  which  exit  in  their  turn.  We 
now  impose  boundary  conditions  at  the  inflow  and  outlfow.  These 
boundary  conditions  correspond  to  the  one-dimensional  conditions 
developed  in  the  earlier  sections,  plus,  in  addition,  the 
condition  v  =  0  at  the  inflow,  and  vx  =  0  at  the  outflow. 

Figure  10  shows  the  influence  of  the  conditions  at  the 
inflow  and  outflow.  The  wave  does  not  continue  as  a  one¬ 
dimensional  wave  but  becomes  smaller  near  the  inflow  and  outflow 
as  the  correct  values  being  imposed  at  these  boundaries  makes 
itself  felt.  Nonetheless,  the  picture  sharply  resembles  the  pure 
one-dimensional  flow  picture  of  figure  9.  There  is  a  small 
amount  of  non-linear  coupling  with  the  other  variables,  but  it  is 
not  significant  enough  to  show  up  on  the  plots. 

The  sidewall  boundary  conditions  used  in  these  plots  were 
the  non-reflecting  characteristic  boundcey  conditions  of  the  last 
section . 

EFFECT  OF  REFLECTING  BOUNDARY  CONDITIONS 

We  now  consider  the  effects  of  other  sidewall  boundary 
conditions.  Figure  11  and  12  show  the  effects  on  T^  if 


reflecting  boundary  conditions  are  imposed  on  the  right 
sidewall.  We  impose  the  conditions  u  =  u,e  =  e,p  =  p  =0 

00  OO  00 

and  Vy  =  0 .  In  terms  of  the  variable  T  ^  ,  this  imposes  the 
condition  T2  +  T3  =  const,  which  ought  cause  a  reflected  wave  in 
T ^  as  the  T2  wave  exits.  Figure  10  corresponds  to  figure  8.  One 
can  see  that  all  variations  are  in  the  y-direction  and  just  as  in 
figure  8,  we  had  the  T2  wave  exit  at  iterations  50-80  so  figure 
10  shows  a  corresponding  T3  wave  entering  and  propogating  into 
the  medium. 

Figure  11  shows  the  corresponding  picuture  for  figure  9, 
when  the  characteristic  boundary  conditions  are  applied  at  the 
inflow  and  outflow.  The  motion  shows  the  influence  of  the  inflow 
and  outflow  boundaries,  so  that  at  some  stage  in  the  future  one 
would  expect  convergence  to  the  free  stream  values,  but  this 
process  is  substantially  showed. 

EFFECTS  OF  OVERSPECIFIED  BOUNDARY  CONDITIONS 

Figure  13  shows  the  interaction  of  the  transverse  wave  with 
the  boundary  when  all  physical  variables  are  specified  at  their 
free-stream  values.  Note  the  presence  of  the  jagged  spikes  as 
socn  as  the  wave  reaches  the  boundary.  This  is  strikingly 
similar  to  the  situation  of  figure  (8)  in  the  one-dimensional 
ca lcula  tion  . 

In  Figure  14,  the  situation  is  almost  exactly  the  same,  the 
principle  difference  being  that  the  overspecified  variables  are 
slightly  different  on  the  two  sidewalls.  This  corresponds  to  the 
situation  where  the  free-stream  values  are  not  known  exactly  but 
only  approximately.  Notice  that  in  the  overspecified  case  of 
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figure  14,  much  larger  errors  are  introduced.  These  errors  only 
seem  to  be  present  near  the  boundary  but  in  more  complicated 
geometries  might  well  propagate  back  into  the  region. 

As  an  example  of  more  serious  problems  with 
overspecification  figure  14  shows  what  ahappens  if  an  initial 
disturbance  of  10%  is  made  in  T2  only.  On  the  left  is  shown  the 
wave  beginning  to  exit,  uneventfully  in  the  case  of  our 
recommended  boundary  conditions.  On  the  right  is  shown  the  large 
errors  occurring  when  the  boundary  is  overspecified.  Note  the 
presence  of  the  large  spikes,  similar  to  the  one-dimensional 
setting  of  figure  7. 

We  conclude  with  a  stiking  example  of  how  one-dimensional 
the  flow  can  be  in  the  entropy  variable  . 

Figure  16  shows  the  effect  of  square  wave  deflection  of  10% 
imposed  uniformly  in  the  x-direction.  Note  how  there  is 
absolutely  no  movement  in  the  y-direction  whatsoever,  although  as 
the  wave  exits,  there  is  predictably  overshoot  visible  at  the  end 


of  the  wave. 


1 


FIGURE  11 

With  reflecting  boundary  conditions  at  sidewall. 


FIGURE  12 


ITERATION  NUMBER  =  ao 


ITERATION  NUMBER 


FIGURE  13 

A  left  running  wave  exits  through  an 
overspecified  boundary-compare  with  figure  8 
(overspecified  1-D)  and  figure  10. 


ITERATION  NUMBER  =  oo  ITERATION  NUMBER 


FIGURE  14 

Another  example  of  overspecification  introduc 
large  spikes-compare  with  figure  13. 


ITERATION  NUMBER  =  10 


iteration  number 


FIGURE  15 

Showing  a  uniformly  wrong  disturbance  in  To  exiting 

a)  on  the  left,  with  characteristic  B.C.'s. 

b)  on  the  righrt,  with  overspecified  B.C.s. 
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FIGURE  16 

Showing  the  surprising  lack  of  transve.  ,e 
movement  in  the  eutropv  variable  T-. . 


CONCLUSIONS 


The  results  of  these  numerical  experiments  show  that  there 
are  four  distinct  artificial  non-reflecting  boundaries  which 
contain  free-stream  information.  These  are  the  inflow,  outflow, 
top-sidewall  and  bottom-sidewall.  Each  of  these  boundaries  must 
be  treated  differently,  or  errors  are  introduced  into  the 
system.  The  worst  errors  are  those  associated  with  over¬ 
specification  and  reflection. 

Other  more  complicated  phenomena  are  present  when  different 
types  of  initial  disturbances  are  present.  All  confirm  the  basic 
lessons  of  this  section,  although  introducing  some  complications 
of  their  own.  At  the  moment,  these  other  results  are  not  fully 
understood  and  call  for  more  experimentation. 

RECOMMENDATIONS 

This  report  raises  many  interesting  questions.  The  most 
obvious  is  that  we  have  investigated  flow  either  perpendicular 
to,  or  parallel  to  an  artificial  boundary,  and  the  boundary 
conditions  required.  However,  in  many  applications,  to  conserve 
the  number  of  grid  points,  the  boundaries  can  be  at  other  angles 
to  the  expected  free-stream  flow.  It  would  be  desirable  to 
deduce  boundary  conditions  where  the  flow  is  at  an  angle  a  to  the 
artificial  boundary,  either  entering  or  exiting  the  region.  It 
is  clear  that  the  conditions,  including,  as  they  would  have  to, 
all  conditions  discussed  here,  would  be  quite  complicated. 

Possible  improvements  to  be  made  in  the  existing  boundary 
conditions  involve  the  improvement  at  the  corners,  There  is  some 


occur 


evidence  of  overspecification  in  the  larger  errors  that 
there . 

It  would  be  interesting  to  compare  these  boundary  conditions 
with  more  complicated  problems  and  possible  oscillatory 


situations . 
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